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Estimation of Purcell factor from mode-splitting spectra in an optical microcavity 
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Washington University, St. Louis, Missouri 63130, USA 

We investigate scattering process in an ultra-high- Q optical microcavity coupled to subwavelength 
scatterers by introducing splitting quality Q sp , a dimensionless parameter defined as the ratio of the 
scatterer-induced mode splitting to the total loss of the coupled system. A simple relation is intro- 
duced to directly estimate the Purcell factor from single-shot measurement of transmission spectrum 
of scatterer-coupled cavity. Experiments with polystyrene (PS) and gold (Au) nanoparticles, Erbium 
ions and Influenza A virions show that Purcell-factor-enhanced preferential tunneling of scattering 
into the cavity mode takes place regardless of the scatterer type. Experimentally determined highest 
<5 ap for single PS and Au nanoparticles are 9.4 and 16.19 corresponding to Purcell factors with lower 
bounds of 353 and 1049, respectively. The highest observed Q Bp was 31.2 for an ensemble of Au 
particles. These values are the highest Q sp and Purcell factors reported up to date. 

PACS numbers: 78.67.Bf, 42.60.Da, 42.81. Qb, 42.65.Es 
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There has been a growing interest in ultra-high- Q op- 
tical whispering gallery mode (WGM) microcavities as 
they provide a suitable platform to study both funda- 
mental physical phenomena (e.g., opto-mechanics, light- 
matter interactions) and to fabricate practical devices 
(e.g., single molecule and nanoparticle detection, opti- 
cal sensors, etc.) [lj. Observed, for the first time, by 
Il'chenko et al. [2| and later addressed comprehensively 
by Weiss et al. [3j, mode-splitting in WGM microcavi- 
ties have developed into a useful tool to investigate light- 
matter interaction [4l4icj. 

Mode splitting occurs due to the lifting of the degen- 
eracy of the clockwise (CW) and the counter-clockwise 
(CCW) propagating WGMs as a result of their coupling 
via scattering centers intrinsically present (i.e., due to 
material inhomogeneity, structural defects, dust contam- 
inations) or intentionally introduced (e.g., particle depo- 
sitions, fiber tips, etc.) into the microcavity mode vol- 
ume. Superposition of CW and CCW modes form two 
orthogonal standing wave modes (SWMs), one of which 
experiences a frequency shift and extra linewidth broad- 
ening with respect to the other. This leads to a doublet 
in the transmission spectrum. 

Previously, we used mode-splitting in passive and ac- 
tive microcavities to measure individual nanoparticles 
[H, EH ■ It was also used to demonstrate Purcell- factor- 
enhanced scattering by considering a collective scatter- 
ing rate of an ensemble of silicon nanocrystals embedded 
randomly into a microcavity Q. These results were in- 
teresting; however they do not provide insight into the 
quantification of mode splitting and the Purcell-factor- 
enhanced scattering at the level of a single scatterer. 

Purcell effect takes place due to light scattering regard- 
less of whether the scatterer is quantum (modeled as dis- 
crete energy levels) or classical (modeled as a dipole). It 
is an important figure of merit for quantifying the ability 
of a cavity to couple to an emitter or a scatterer. Thus, 
it is important to develop reliable and easily accessible 



methods to measure Purcell factor accurately. In this 
Letter, we report the direct estimation of Purcell factor 
from mode splitting in the transmission spectrum with a 
single-shot measurement, and clarify the role of individ- 
ual scatterer on the quality of mode-splitting and hence 
on Purcell-enhanced Rayleigh scattering. 

The proposed method for direct estimation of Purcell 
factor and demonstrated high Purcell-enhanced scatter- 
ing efficiencies at single scatterer resolution, using PS and 
Au nanoparticles, and Influenza A (InfA) virions, will 
play an important role in detecting and studying clas- 
sical and quantum scatterers independent of their type 
and internal structures. For example, quantum scatter- 
ers such as atoms and molecules can be detected by mode 
splitting spectra and Purcell-enhanced scattering in a mi- 
crocavity with a resonance largely detuned from their op- 
tical transitions such that they are not optically excited. 
This will provide an alternative technique to fluorescence, 
absorption and ionization based atom detection and will 
eliminate the need for near-resonant lasers. 

The amount of mode-splitting and the difference in the 
linewidths of the individual resonances in the doublet 
are given as 2g — —af 2 (r)oj/V and 2T = —2gauj 3 /3irc 3 
where the polarizability a is defined as a — 4irR 3 (n 2 — 
1)/ (n 2 + 2) for a single particle of radius R and refractive 
index n, V denotes the mode volume, f(r) the normal- 
ized mode distribution, c the speed of light, and w the 
angular frequency of the initial resonant WGM 0, H[. 
The coupling strength of the initially degenerate WGMs 
can be described with the mode- splitting quality Q sp as 

Q - 2|g| (i) 
Wsp r + ( w /Q) {) 

where w / Q is the initial resonance linewidth. To resolve 
the mode-splitting in the transmission spectra, Q sp > 1 
should be satisfied, i.e., 2\g\ should exceed the mean of 
intrinsic and scatterer induced decay rates of the system: 



2\g\ >r+(«/Q). 



(2) 
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FIG. 1: Splitting quality versus particle size (R) for PS (re- 
fractive index n = 1.59) particles for initial cavity quality 
factor Q equals to 10 6 and 10 8 . Vertical red lines denotes the 
size range of particles used in the experiments. Mode split- 
ting is resolved if Q sp lies above the horizontal dashed line 
(Q Bp > 1). The points labeled with * denote the Q sp values 
obtained in single-particle experiments with particles of mean 
radii R = lOOnm and R — 135nm. 



If the scattering losses dominate the cavity losses, i.e., 
r > (w/Q), Eq. flTJ simplifies to $ sp = 2\g\/T. Then 
quality of mode-splitting can be related to the nanopar- 
ticle through particle polarizability as 



Q\r 
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Thus, in this regime mode-splitting quality is affected 
neither by the Q of the cavity modes nor by f(r). It 
is determined by the resonance wavelength of the initial 
WGM and the particle polarizability. 

On the other hand, if the cavity losses dominates, i.e., 
r < (u/Q), Eq. (HJ simplifies to 



2\g\Q _af 2 (r)Q aQ 
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where the upper bound is obtained by setting f 2 (r) = 1. 
Thus, in this regime Q sp increases linearly with particle 
polarizability and the quality factor of the WGM. More- 
over, the mode field distribution and volume strongly af- 
fect the splitting-quality in this regime. Since one does 
not have precise control on the location of the nanoparti- 
cle within the microcavity mode volume, f(r) is unknown 
and it plays a limiting role for <3 sp in this regime. The 
expressions obtained for splitting quality for scattering 
losses dominant and for cavity losses dominant regimes 
in Eqs. (|3]) and (|4|), respectively, also establish the rela- 
tionship for the lower bound of Purcell factor T: 



(5) 



where T = (3/47r 2 )A 3 (<2/V). The second part of the in- 
equality is obtained using l/Q sp = 1/Q^ P + l/Qsp an d 
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FIG. 2: Total quality factor Q T versus mode splitting quality 
Q sp for various values of PS particle size 7?. Red * points 
correspond to randomly chosen (Q T ,R) pairs in the region 
10 6 < Q < 5 x 10 7 and 50 nm < R < 100 nm. Linear 
boundaries are obtained by varying Q in 10 6 < Q T < 10 s for 
fixed R, and inverted parabola-like boundaries are obtained 
by decreasing R down to 10 nm for fixed Q. 



to estimate the lower bound of Purcell factor in an ex- 
periment by using the <3 sp obtained from transmission 
spectrum of single-particle coupled microcavity. 

In Fig. [TJ we depict the dependence of Q sp on the par- 
ticle size and the quality-factor of the microcavity. For 
particles larger than a critical size, we have <3 sp ~ Q l sp 
suggesting that the additional decay channel formed by 
the particle dominates the cavity losses. On the other 
hand, for particles smaller than a critical size, we see 
Qsp ~ Qlp which suggests that small particles do not in- 
duce significant loss when compared to the cavity losses. 
In between these two critical particle sizes, there is a 
region where the splitting-quality can be explained nei- 
ther by Q 1 nor by Q^ p ; hence the exact expression of Q sp 
given in Eq.([rj should be used. Furthermore, we see that 
for a given Q, there is a critical lower and critical upper 
bound for the particle size beyond which Q sp becomes 
less than one. While the upper bound is determined by 
particle size, the lower bound is strongly affected by Q. 
Therefore, by choosing a mode with higher Q, mode- 
splitting induced by particles with radii of a few nanome- 
ters can be resolved, subsequently such particles can be 
detected and measured. It is worth noting that Q sp for 
a fixed cavity Q- factor is upper bounded by C^ sp = <^ p , 
i.e., uj/Q =T implying Q sp = Q* p /2 = Q^/2, and the 
particle size which pushes <5 sp to its maximum value be- 
comes smaller with increasing cavity Q-factor. One can 
easily show that the particle maximizing <5 sp should have 
polarizability satisfying a = 3n~ 2 f (r)^ 1 \ 3 ^/T / 4. 

Here we discuss how the mode-splitting quality <3 sp 
depends on the total quality factor Q T of a scatterer- 
coupled microcavity (FigfJ]). Q T takes into account both 
the cavity losses and particle related losses, i.e.,Q x = 
w[r+ (uj/Q)}- 1 . Then Q sp = Q T (2g/uj) implying that 
<3 sp increases linearly with Q T with a slope of 2g/u if 
R (or a) is kept constant, i.e., g and T are fixed too. 
When Q is very small, process is dominated by cavity 
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FIG. 3: Schematics of a tapered fibre interfaced to a silica mi- 
crotoroid with deposited nanoparticles, and typical transmis- 
sion spectra captured with a photo detector for PS particles 
of radii R = lOOnm (top) and R = 135nm (bottom). 



For multi-scatterer case, SWMs are affected by all par- 
ticles in varying degrees and they are spectrally-shifted 
from the original WGM resonance. We denote (cOj^Jj) 
with j = 0, 1, 2 as the angular resonance frequencies and 
linewidths of the initial (pre-scatterer) WGM (j = 0), 
and low (j — 1) and high (j — 2) frequency SWMs cre- 
ated after the deposition of N scatterers. Then, we have 
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losses. Since scatterer-related-losses are constant (i.e., 
R is fixed), increasing Q decreases cavity losses. When 
cavity and scatterer losses are equal (r = uj / Q) , we have 
Q T = uj/2T and <5 sp = <3 sp /2 = <?/T. If Q is increased 
beyond this value, cavity losses becomes smaller than the 
scatterer losses. In the limit of large Q, we have Q T = 
uj/Y leading to <5 sp = Q\ p = 2g/T. Thus two sections can 
be defined: (LI) Cavity losses dominant section (uj/Q > 
r) which is the part of the linear curve within < Q T < 
uj/2T and < Q sp < Q\ p /2. (L2) Scattering dominant 
section (uj/Q < V) is the part bounded with Q T > u/2T 
and Q* p /2 < Q sp < Qi p . 

Next we investigate what happens if Q is kept constant 
at Q and the particle size R is decreased from its ini- 
tial value of i?o (FigHJ). If the process is initially in the 
region oj/Q > T Ro , decreasing particle size will signifi- 
cantly decrease V due to its V cx R 6 dependence which 
makes it much lower than oj/ Q . Due to g oc R 3 depen- 
dence, g also decreases but at a slower pace. Then, in 
this case the process continues to be dominated by cav- 
ity losses, and Q sp decreases with decreasing R. If the 
process is initially in the region uj/Q < ^r , the evo- 
lution is strongly affected by g and T rather than uj/Q. 
Thus, the slower pace of decrease in g compared to that 
in r allows an increase in Q sp . This continues until the 
point of ijj/Q = Tji and <5 sp = g/T. Further decrease of 
R shifts the process from scatterer-dominated regime to 
cavity- loss-dominated regime in which the effect of oj/Q 
is much stronger. Thus, after stepping into this regime, 
the process will stay in this regime with decreasing R 
until <3 sp decreases from its peak value to zero. 

Figure [3] shows transmission spectra obtained for single 
PS particles of R = lOOnm and 135nm coupled to a mi- 
crotoroid [f| . For these spectra, we calculate <5 sp = 8 and 
Q sp = 9.4 from Eq. (P). Then Eq. © yields the corre- 
sponding Purcell factors as F > 256 and T > 353. This is 
the best Purcell factor obtained for a single PS nanoparti- 
cle. The * points in FigJTJdenote experimentally obtained 
<5 sp for different depositions of nanoparticles whose sizes 
were estimated using a = (3A 3 /87r 2 )(l/(3^ p ) 6]. Differ- 
ent Q sp values for the same estimated size are due to 
different f(r) at each deposition. 



where 4>i denotes the spatial distance of i-th particle from 
the anti-node of a SWM. From Eq. ([6]), we find the 
amount of mode splitting as 



N 



X = U2 - ui = ^ 2ffj cos(2(fo 



Moreover, Eq. ([7]) yields the linewidth relation as 



(8) 
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7i + 72 = 270 + ^2 2T t 
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To resolve the split modes after the deposition of N scat- 
terers, |x| > (71 +72V2 should be satisfied. Thus, mode 
splitting quality for TV-scatterer then becomes 



j(N) 



2|X| 



Eili 9% cos(20i) 



71+72 (wo/Q) + EL r < 



(10) 



which has a form similar to Eq. (fT]). Assuming that the 
particles have the same polarizability and using triangle 
inequality and < |cos(20j)| < 1, we find that split- 
ting quality in the limit of scatterer dominant regime 
becomes Q^ 1 < Ql p where Q 1 is given in Eq. ([3]). 
If we further set / 2 (?'i) = / 2 (r) and consider the cavity 
losses dominant regime, mode splitting quality becomes 
Q ( S p )n = af 2 {r)QN/V < aQN/V where in the last part 
we assumed maximum overlap between the WGM and 
each of the scatterers, i.e.,/ 2 (r) = 1. Assuming that all 
particles have the same a, and defining g = g 1 f 2 (ri) and 
r = T'f 2 (n) with g' = -auj/2V and V = a 2 w 4 / 'GirV c 3 , 
it is easy to see that 
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where we used ff{r) — f 2 (r) = 1. For N = 1, Eq. 
(fTTj) reproduces the results for the single scatterer case. 
Following the procedure carried out for a single scatterer, 
we find the Purcell factor for multi-scatterer case as 
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(12) 



4 




! 300 



&■ 200 
S 

S" loo 
£ o 



X 15 
£ 10 

s 

1 5 

5 o 



10 15 20 25 30 35 40 45 



5 10 15 20 25 30 35 40 45 
Time (s) 



FIG. 4: A series of transmission spectra obtained for consec- 
utive depositions of PS nanoparticles (R = 135nm) into the 
mode volume of the resonator and the corresponding amount 
of mode splitting, additional linewidth broadening and the 
mode splitting quality. A total of five nanoparticle deposition 
is depicted. Each frequency jump in the spectra corresponds 
to a single nanoparticle deposition event. 



from which one can deduce that the minimum attain- 
able Purcell factor decreases with increasing number of 
scatterers due to the increasingly dominating scatterer 
induced dissipation. Thus lower bound of the Purcell 
factor can be estimated from the experimentally obtained 
transmission spectra and using the calculated Q^p^ in Eq. 
(|12[) provided that the number of particles TV is known. 

Figure 2] shows the evolution of transmission spec- 
tra and calculated mode splitting, additional linewidth 
broadening and mode splitting quality as nanoparticles 
are deposited into the resonator mode volume. Each dis- 
crete jump corresponds to deposition of a single nanopar- 
ticle. With each particle deposition, mode splitting is 
either enhanced or reduced depending on the location 
of the particle on the mode volume. Measured values 
of QW = 9.2, qP) = U.I, Qg) = 14.0, G$ - 7.3 
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and = 6.8, respectively, corresponds to J 7 ^ > 338, 
F^ > 397, F^ > 261, > 53, and J*® > 36. Esti- 
mated Purcell factor for two- particle event, which is also 
the highest Purcell factor ever reported d, [l2| , implies 
that at least 99.75% of the scattered light is captured by 
the doubly degenerate cavity modes [8| . 

We investigated the effect of scatterer properties with 
experiments using Au nanoparticles, Erbium (Er +3 ) ions 
and InfA virions. For a single Au nanoparticle (R = 
lOOnm), we measured Q sp = 16.2 corresponding to 
F^ 1 ' > 1049. With subsequent depositions, <5 sp changed 
as = 20.16, 0$ = 24.64 and Q$ = 20.88 with 



the estimated F {2 ^ > 812, F { ^ > 809, and F^ > 435. 
For an ensemble of N < 20, the highest measured Q sp 
was 31.2 with the resonance modes shifted from each 
other by 2\g\ = 364.9MHz with a mean linewidth of 
11.7MHz. This corresponds to F > 194.7 implying a cap- 
ture efficiency higher than 99.5%. The splitting quality 
of <5 sp = 31.2 is the highest reported up to date [8j. Next, 
we performed experiments using Er-doped silica micro- 
toroids [Til] - Gain provided by optically pumping Er +3 
ions enables detection of small intrinsic mode splitting, 
which could not be observed before optical pumping. 
Such small mode splittings yielded <5 sp ~ 3.5. Finally, we 
deposited InfA virions one-by-one on a microtoroid. Esti- 
mated mode splitting qualities are = 2.5, = 5.9, 
Qi 3 p> = 9.8, q£> = 9.2 and Q$ = 12.6, respectively 
corresponding to F {1 ~> > 25, J" (2) > 69, J" (3) > 128, 
F^ > 84, and F {5 ~> > 127. These results demonstrate 
that preferential backscattering into the resonator takes 
place regardless of the scatterer type, and that it can 
be quantified from the mode splitting spectrum with a 
single-shot measurement. 

In conclusion, we have investigated scattering induced 
mode splitting in ultra-high- Q microcavities, and intro- 
duced a formalism to quantify the quality of mode split- 
ting for single and multiple-scatterers as the particles 
enter the mode volume of the resonator. This formal- 
ism enables estimating the lower bound for Purcell fac- 
tor and Rayleigh scattering capture efficiency from the 
transmission spectra. This will help understanding the 
Purcell-factor-enhanced scattering process at the single 
particle level, and boosting the sensitivity of detecting 
individual nanoparticles entering to or embedded into 
the mode volume of a resonator. Finally, this method 
can be further extended for a wide class of microcavi- 
ties with embedded, adsorbed or deposited emitters and 
scatterers (e.g., nanocrystals, plasmonic particles, atoms, 
molecular gases) independent of their internal structures 
using collective or individual enhancement effects, thus 
providing a non-invasive detection scheme. 
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